Tomography and state reconstruction with superconducting single-photon detectors 
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We investigate the performance of a single-element superconducting single-photon detector 
(SSPD) for quantum state reconstruction. We perform quantum state reconstruction, using the 
measured photon counting behavior of the detector. Standard quantum state reconstruction as- 
sumes a linear response; this simple model fails for SSPDs, which are known to show a non-linear 
response intrinsic to the detection mechanism. We quantify the photon counting behaviour of the 
SSPD by a sparsity-based detector tomography technique and use this to perform quantum state 
reconstruction of both thermal and coherent states. We find that the nonlinearities inherent in the 
detection process enhance the ability of the detector to do state reconstruction compared to a linear 
detector with similar efficiency for detecting single photons. 



I. INTRODUCTION 

We investigate the photon counting abilities of Super- 
conducting Single-Photon Detectors (SSPDs) (lijjj . mo- 
tivated by the recent u psu rge in the use of such detec- 
tors in quan tum optics [19 . |20| and quantum cryptogra- 
phy pll |22|| experiments. SSPDs are fast [23], spectrally 
broadband |24|-l26j single-photon sensitive detector with 
low noise [27| . These detectors consist of an ultrathin me- 
andering strip of a superconductor with low Cooper pair 
density, typically NbN. When biased close to the critical 
current, the absorption of a photon produces a transition 
from the superconducting to the normal state, resulting 
in the creation of a resistive area and the appearance of 
a voltage pulse in the external read-out circuit. 

It was previously shown 0, [3 [H, [H, [H| that de- 
pending on the bias current through the superconductor, 
the detector has multiphoton regimes, where the energy 
from several photons is required to break the supercon- 
ductivity. These multiphoton detection events - which 
depend on several photons being absorbed close together 
- increase the probability of the detector clicking in a 

nontrivial way Em? ]. 

Quantum state reconstruction involves finding the pho- 
ton number distribution of an unknown quantum state of 
light from the response of some detector. This task is of 
fundamental importance for any quantum optics or quan- 
tum communication experiment, as the final step in such 
an experiment is always the measurement of the photon 
occupation number in a particular detection mode. Sur- 
prisingly, the reconstruction of a radiation field can be 
one with a single detector that has only an on/off output 
This is possible because measuring the count rate at 
different settings of the detector (e.g. at different efficien- 
cies) produces a response curve characteristic for each 
photon number. The reconstructed distribution of pho- 
ton numbers for an unknown state is then given by the 
linear combination of response curves that best describes 
the measured count rates of that state d, H, [l2rlT7| . By 
taking into account the finite efficiency of the detector, 



it is possible to reconstruct the state at the input of the 
detector rather than the statistics of the absorbed pho- 
tons. 

The results in this paper are divided into three parts. 
In Section III, we perform a modified version of detector 
tomography [7] on the SSPD. This tomography quanti- 
fies the complex behaviour of the device, enabling the 
use of the SSPD in situations where responses to sev- 
eral different photon numbers are important. By using 
a technique which has minimal assumptions, we over- 
come the problem that the understanding of the work- 
ing of the detector is still incomplete [30] , harnessing the 
SSPD for quantitative multiphoton applications. Our to- 
mographic technique is based on sparsity. The advantage 
of a sparsity-based technique is its robustness [29|. We 
show explictly how to apply this tomographic technique 
to an SSPD. 

Next, in Section IV, we perform quantum state recon- 
struction with the SSPD. We reconstruct states with av- 
erage photon numbers up to (n) = 11.4, thereby showing 
that the tomographic process was succesful, and demon- 
strating that it is possible to use a nonlinear device for 
quantum state reconstruction. We reconstruct both co- 
herent and thermal states. 

Lastly, in Section V, we evaluate the effect of these 
nonlinearities on the quality of the state reconstruction 
process. By evaluating the Cramer-Rao bound - the the- 
oretical limit on the amount of information that may be 
extracted from a measurement - we establish that the in- 
trinsic nonlinearities of an SSPD are benificial for quan- 
tum state reconstruction. 



II. THEORY: QUANTUM STATE 
RECONSTRUCTION 

The goal of Quantum State Reconstruction is to recon- 
struct the diagonal elements p nn = diag(p) of the input 
state of the light. From observations at different config- 
urations of the detector, it is possible to reconstruct the 
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state because each photon number gives a particular re- 
sponse on the detector that is characteristic of that pho- 
ton number. The task is to find from the set of response 
curves the linear combination of photon numbers that 
best describes the measured count rate of some unknown 
state. The response of the detector at each setting v of 
the tuning parameter is described by a Positive Operator- 
Valued Measure (POVM) element n„ = £„IL„ \n)(n\, 
and the detector responds to the state with a probability 
R v = Tr(pU v ). The POVM contains a full quantitative 
description of the measurement process. 

Due to the shot noise associated with the discreteness 
of the photon counting process, the problem of solving 
this set of equations simulataneously is inevitably statis- 
tical in nature, since the equations will not be analyticaly 
invertible. This problem can be solved by a maximum 
likelihood (ML) technique, usin g th e Expectation Max- 
imization (EM) algorithm [f| |R1T^[l7jj to find the best 
solution, while respecting the normalization of the state. 
A derivation is given in [l5[. The i-th iteration of this 
algorithm is given by: 

Ji+D = Ji) >p LLvn n » m 

where pW is the state at the i-th iteration, N is the total 
number of experimental preparations, R v is the measured 
click probability at the v-th experimental configuration 
and Pi/(p^) is the calculated click probability at the z^-th 
experimental configuration. It is known that this algo- 
rithm converges to the ML solution, for which the stan- 
dard errors are given by the Cramer- Rao bound [17]. It is 
also known that this algorithm can take many iterations 
to converge. Following earlier work pj M. fl^ - [l7| . we take 
our number of iterations to be 10 6 . 



III. EXPERIMENTAL SETUP 

The SSPD used in this experiment is a commercial 
NbN meander produced by Scontel. The width of the 
wire is 100 nm, and the distance between the wires is 
150 nm. The size of the active area is 10 pm by 10 pm. 
The device was cooled in a bath cryostat to a temper- 
ature of 1.7 K. The measured overall system quantum 
efficiency for the one-photon Fock state was 2.8% at a 
bias current of 13.3 pA (corresponding to 7&/I c « 0.9) 
and a wavelength of A = 1500 nm. 

For our detector tomography procedure, we illuminate 
the device with a series of coherent states varying from 
130 fW to 108 nW (0.05 to 4.1*10 4 photons/pulse). The 
low powers were achieved with a computer-controlled 
variable attenuator, whose linearity to -60 dB was verified 
independently. From the measured response to coherent 
states, we reconstruct the POVM using the method de- 
scribed below. The coherent states were generated by 
a Fianium supercontinuum pulsed laser. The repetition 



rate of this laser was 20 MHz, the specified pulse width < 
7 ps. The light was filtered to have a center wavelength 
Ao = 1500 nm and a spectral width AA = 12 nm. The 
observed POVM was then used to reconstruct coherent 
and thermal states. We verified independently that the 
output from our supercontinuum laser is indeed a coher- 
ent state. We measure g^ 2 '(0) with a coincidence circuit 
and obtain g {2) {0) = 0.97 ± 0.02. 

We generate pseudothermal states by the standard 
technique of a rotating ground glass plate [32l |. which 
was illuminated with the coherent states described above. 
The exponential probability distribution of the intensity 
of the resulting speckles creates photon statistics that 
are equivalent to thermal light when averaged over many 
realizations of the angle setting of the plate. 

Unfortunately, after the reconstruction of the coherent 
states, the alignment of the detector in the cryostat was 
degraded. We therefore recharacterized the device in its 
new configuration with a set of coherent states before 
performing the reconstruction of the thermal states. The 
degradation manifests itself as an increased dark count 
probability, which was 0.01 / pulse at h/I c ~ 0.9. 



IV. SPARSITY-BASED TOMOGRAPHY 

We start by measuring the detector response curves, 
i.e. the detection probablity versus detector bias cur- 
rent, for a set of coherent states. From these, we deduce 
the more fundamental response curves for Fock states. 
Fig. 1 shows the resulting set of inferred detector re- 
sponse curves, i.e. the probability of the detector to re- 
spond to a certain Fock state. The detector tomography 
is per formed by a method based on the one described 
in [29]. The essential assumption is one of sparsity: we 
describe the detector by as few physical parameters as 
possible, while not restricting the possible range of be- 
haviours that our model describes. More specifically, we 
describe our detector by a combination of linear attenua- 
tion (given by linear losses in the detection process such 
as the finite absorption of the NbN layer) followed by a 
nonlinear photodetection process inside the NbN layer. 
The reason for including a linear efficiency seperately is 
that it significantly reduces the number of parameters 
required to model the detector, making the tomography 
more robust. 

A second reason is that nonunity linear efficiency in- 
troduces correlations between the variousn„„ at one bias 
current. The reason for this is that at nonunity efficiency, 
the n photons necessary for an n-photon nonlinear pro- 
cess could have come from any N > n number of incident 
photons. By explicitly including this effect, we make sure 
that our reconstructed POVM is compatible with this 
process. 

This description in terms of a linear absorber and a 
nonlinear process, which we showed to be applicable to 
the NbN nanodetector (3lj is applicable to the SSPD as 
well. We can therefore write for the click probability: 



3 



k—oo 
k=0 



k\ 



(2) 



where r\ is the linear efficiency, (n) is the mean photon 
number, R u is the count rate at a given bias current and 
the p v k are the POVM elements prior to the inclusion of 
the finite linear efficiency. These elements now only in- 
clude the nonlinear effects of the detector, i.e. the photon 
number threshold regime that the detector is in, which 
depends on the bias current. 

After each fit, where the fits at different currents are 
completely independent, we reinclude the r\ into the p u k 
to produce the POVM element Ii u k by the following pro- 
cedure: first, we fix the number n mr at which we are go- 
ing to truncate the Hilbert space for the reconstruction. 
Then, we construct a vector of length n mr , where the first 
5 elements are and the other elements are equal to 
1. Finally, we multiply this vector by a Bernouilli trans- 
formation [? ] Lkk' = ( fei ^ V k (1 — r)) k ~ k , absorbing 

the linear losses into the POVM. We perform state recon- 
struction with the POVM consisting of all II„fc obtained 
at different currents. 

For this experiment, we are not interested in separating 
the linear and nonlinear effects, but rather in finding a 
description of the entire detector. Therefore, we truncate 
the sum at n max = 4 for all currents. This is equivalent 
to assuming that the detector is governed only by linear 
effects at sufficiently high photon numbers. This choice 
is motivated by the fact that we do not enter the three- 
photon regime in the current range over which we operate 
our detector and is justified by the good fits obtained with 
this model. 

For the analysis, we grid our measured count rates by 
linear interpolation, producing 165 current settings, from 
5 /iA to 13.25 fiA, which is the current range over which 
we could measure count rates at enough powers to create 
a good fit. This is also the current range over which we 
perform the state reconstruction. 



V. QUANTUM STATE RECONSTRUCTION 

Fig 2. shows a representative sample of the recon- 
structed coherent and thermal states. We reconstruct 
a series of coherent and thermal states, using the al- 
gorithm given by Eq. (1), iterated 10 6 times. For the 
quality of our reconstruction we use the fidelity, defined 
as F — X)"=o° \/ PnnPnn, where p represents the density 
matrix of the coherent state corresponding to the average 
number of photons found in the reconstruction. 

Fig. 3 shows the fidelity of the state reconstruction, as 
a function of mean photon number (n). We observe that 
the quality of the reconstruction degrades as the average 
number of photons increases. This can be understood 
from Fig. 1: as the number of photons increases, the 
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Figure 1: Response curves inferred from detector tomography 
as a function of bias current through the device. On the x axis 
is the bias current, on the y axis is the probability that the 
detector responds to a particular number of photons (Fock 
state). The black line indicates photons, the arrow indicates 
the direction of increasing photon number. Note that we have 
shown only the first incident 15 photon numbers for clarity. 
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Figure 2: Six typical reconstructed states. a)-d) show coher- 
ent states, e) and f) are thermal states. The red bars show the 
closest coherent (for a-d) or thermal (e & f) state, the black 
bars show the reconstructed state, i.e. the result of eq. 1 after 
10 6 iterations. The fidelities and mean photon numbers are 
indicated for each reconstructed state 



4 




■ coherent 

theory coherent 

O thermal 

- ■ theory thermal 



O 

o 



1.00 
0.98 
0.96 
0) 0.94 
^ 0.92 
0.90 
0.88 

2 4 6 8 10 12 

mean photon number <n> 

Figure 3: Fidelity of the reconstructed states as function of 
mean photon number. The solid black squares indicate recon- 
structed coherent states, the open red circles indicate recon- 
structed thermal states. The lines indicate the results for the 
expected fidelities from simulations. 



curves lie closer together, making it more difficult to dis- 
tinguish the contributions from various photon numbers. 

The theoretical curves in Fig. 3 were generated by sim- 
ulating the experiment. Each experiment was simulated 
30 times, to obtain a reasonable estimate of the expected 
fidelity. The simulations were performed by calculating 
expected count rates from the POVM and a given state. 
For each calculated count rate curve we assumed a con- 
stant relative error. We made the approximation that 
these errors are uniformly distrubuted within the inter- 
val [R-AR, R + AR}. 

We compared the mean square error x 2 for the re- 
constructed coherent and thermal states with the the- 
oretically expected count rates for coherent and thermal 
statistics. From this analysis we conclude that we can 
successfully distinguish between coherent and thermal 
states. At larger values of (n) the value of \ 2 becomes 
large, indicating that the state reconstruction becomes 
inaccurate and loses its ability to correctly predict the 
quantum state. This happens at (n) w 9 and (n) « 15 
for the thermal and coherent states respectively. 

By comparison with the measurements, we find that 
the relative error AR/R is 2% for the coherent states, 
and 6% for the thermal states. These numbers are justi- 
fied by the observed deviations between count rates ex- 
pected from the reconstructed states and the measured 
count rates. We attribute this error to the uncertainty 
in setting the bias current through the device, where 
AR/R = 2% corresponds to 40 nA of uncertainty in the 
bias current. We attribute the higher uncertainty for the 
thermal states to residual variations in the input inten- 
sity caused by the rotating ground glass plate. 
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Figure 4: Theoretical limitations on state reconstruction, as 
calculated from the Cramer-Rao bound, for the same coherent 
state ((n) — 2.5) measured by two different detectors. This 
figure shows the limits on the ability of the detector to find 
the components p„ n of this state. We show the relative error, 
i.e. 5 p nn = <J Pnrl /pnn- The black points indicate the errors 
expected for our SSPD, the red ones are for an APD of linear 
efficiency equal to our SSPD, where the linear efficiency is the 
tuning parameter. The errors are calculated assuming 6 * 10 s 
measurements divided over 100 settings of the tuning param- 
eter, corresponding to 30 s experimentation time. This figure 
shows that the nonlinear effects in the SSPD are beneficial for 
state reconstruction. 



VI. NONLINEARITY-ENHANCED QSR 

In Fig. 4, we show the effect of the Poissonian (shot- 
noise) errors on our reconstruction, calculated from the 
Cramer-Rao bound [§1, In this figure, we compare 
the expected error bars for state reconstruction of our 
SSPD with those of an APD of linear efficiency equal to 
the SSPD. The shot-noise error is the fundamental lower 
limit on the error and is fixed for a given state. There- 
fore, it is relevant to investigate the state reconstruction 
abilities of a detector. 

Fig. 4 shows that the errors in the reconstruction of the 
SSPD are 50% lower than those of an APD with equiva- 
lent efficiency. We attribute this to the nonlinear effects 
which give our detector higher efficiency at the multipho- 
ton level [i[ . The physical reason for this nonlinearity is 
that two photons which are absorbed close together on 
the nanowire have an enhanced probability to make the 
detector click. This is in contrast with an avalanche pho- 
todiode, where - as long as each photon overcomes the 
bandgap - there is no mechanism where the photons as- 
sist one another in producing an avalanche. 



VII. DISCUSSION 

In previous state reconstruction experiments, 
avalanche photodiodes (APDs) were used as detec- 
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tors Q , where the tuning parameter was the attenuation 
of an extra attenuator inserted in front of the detector. 
These detectors have a POVM element determined 
only by linear attenuation, combined with single-photon 
threshold behaviour. Banaszek has noted before (l5j 
that state reconstruction is not limited to T\ vn containing 
only linear attenuation. To our knowledge this is the 
first time that state reconstruction has been performed 
with a POVM that contains nonlinear as well as linear 
terms. 

We find that the practical limitations of our experi- 
ment are the accuracy with which we can set the current. 
We note, however, that this is no fundamental limita- 
tion, as there are current meters available which have a 
much higher resolution than the one used in our experi- 
ment. Furthermore, we note that the Cramer-Rao bound 
(CRB) has the usual square-root dependence on the num- 
ber of measurements made, indicating that longer mea- 
surement times will improve the quality of the reconstruc- 
tion as well. 

The EM algorithm is known to converge to a solution 
that saturates the CRB, meaning that we achieve the 
lowest possible variance (i.e. the optimal solution) for 
our reconstruction. Moreover, through propagation of 
errors, the limits on state reconstruction set hard limits 
on the ability of the detector to measure other quantitites 
which are functions of the input state, such as the second- 
order correlation function g^ 2 '(0) [33||. Therefore, state 
reconstruction is a fundamental tool to investigate detec- 
tors operating at the few-photon level, e.g. multi-element 
SSPD detectors [28l l3343q . Since state reconstruction is 
performed at the limit of the amount of information that 
can be extracted from a measurement, investigating the 
state reconstruction abilities of a detector may yield bet- 
ter understanding of its capabilities. Furthermore, since 
QSR probes the limits of the abilities of a detector, we 
propose it as a useful benchmark tool to compare various 
SSPDs. 

In this work, we have used a commercial SSPD, with 
an efficiency of up to 2.8%. However, recently, great 
strides have been taken in making SSPDs more efficient, 
by including cavity structures on the SSPD [37l ? ? ], by 
switching to different materials such as WSi [? ] and 
by incorporation of SSPDs in nanophotonic waveguides 
[? ? ]. Marsili et al report efficiencies as high as 93%. 
The present work opens up the possibility of using such 
high-efficiency devices for quantum state reconstruction. 

It is an open question whether these more efficient de- 
tectors would have beneficial multiphoton nonlinearities 
similar to the detector reported on in the present work. 



We note that this nonlinear effect has also been reported 
in |4| , and since it is due to a physical effect that will also 
be present in more efficient detectors (namely the absorp- 
tion of two photons close together) there is nothing that 
precludes it. Detector tomography of a high-efficiency 
SSPD is necessary to answer this question. 

The present work opens up the possibility of applying 
the ideas of optimal experimental design to the design of 
SSPDs. For quantum tomography of a spin-1/2 system, 
there has been work [39j on how to optimize the POVM 
used in a measurement to yield the optimal tomogra- 
phy result. Similar reasoning may be applied to design 
an SSPD that is particularly suitable for reconstructing 
particular states, or for measuring particular properties 
of states. Such optimization - constrained by what is 
possible with present production techniques - would fo- 
cus on the width of the wire (which governs the nonlin- 
earity), length of wire segments, number of elements in a 
multi-element detector, and the relative size and detec- 
tion efficiency of each element. 



VIII. CONCLUSIONS 

We have shown quantum state reconstruction with an 
SSPD. The SSPD is especially suited for this task because 
of its low noise, fast response time, its intrinsic tuning 
parameter in the form of the bias current through the 
device, and the nonlinearities which enhances the state 
reconstruction abilities. 

Since the fundamental physics of this device is not 
known, we have performed a detector tomography pro- 
cedure in order to find the parameters describing the 
response of this device (POVM). We thus demonstrate 
state reconstruction with an experimentally determined 
POVM. This illustrates the utility of detector tomogra- 
phy- 
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